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Abstract 

We present the first public code for semi-analytical calculation of the y-ray flux astrophysical /-factor from dark matter annihila- 
tion/decay in the Galaxy, including dark matter substructures. The core of the code is the calculation of the line of sight integral of 
the dark matter density squared (for annihilations) or density (for decaying dark matter). The code can be used in three modes: i) 
to draw skymaps from the Galactic smooth component and/or the substructure contributions, ii) to calculate the flux from a specific 
halo (that is not the Galactic halo, e.g. dwarf spheroidal galaxies) or iii) to perform simple statistical operations from a list of al- 
lowed DM profiles for a given object. Extragalactic contributions and other tracers of DM annihilation (e.g. positrons, antiprotons) 
will be included in a second release. 
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^ PROGRAM SUMMARY 

^.Program Title: CLUMPY 

Q Programming language: C/C++ 

Computer: PC and Mac 
CZ2 .Operating system: UNIX(Linux), MacOS X 
, ^ 'RAM: depends on the requested size of skymaps (~40 Mb for a 

|500 x 500 map) 
y—( -Keywords: dark matter, indirect detection, gamma-rays 
Classification: 1.1, 1.9 

External routines/libraries: CERN ROOlQ, Doxygei£] (optional) 
Nature of problem: Calculation of y-ray signal from dark matter an- 
nihilation (resp. decay). This involves a particle physics term and an 
astrophysical one. The focus here is on the latter. 
Solution method: Integration of the DM density squared (resp. den- 
sity) along a line of sight. The code is optimised to deal with the DM 
density peaks encountered along the line of sight (DM substructures). 
"A semi-analytical approach (calibrated on N-body simulations) is used 
,for the spatial and mass distributions of the dark matter substructures 
• i-h >in the Galaxy. 

Restrictions: Some generic dark matter annihilation spectra are pro- 
vided but are not included in the calculation so far as it is assumed that 
the particle physics is independent of the astrophysics of the problem. 
Running time: this is highly dependent of the DM profiles considered, 
the requested precision e and integration angle arjn t : 

• about 60 mn for a 5° x 5° map towards the Galactic centre, with 
a mt = 0.01°, NFW dark matter profiles and e = 10~ 2 ; 

• about 2h for the same set-up towards the anti-centre; 

• 0.1 to 10 DM models per second, depending on integration angle 
and DM profile. 
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1. Introduction 



Despite the several astrophysical evidences pointing at the 
existence of dark matter (flat rotation curves of galaxies, gravi- 
tational lensing, "bullet cluster", etc.), its nature still evades us. 
This search has become one of the major topic in both particle 
physics and astrophysics and is tackled using either direct or in- 
direct detection methods. For the former, the hope is to directly 
witness the interaction of a dark matter particle with a detec- 
tor. The indirect approaches aim at measuring the end products 
of dark matter annihilation/decay (e.g., y-rays, positrons, anti- 
protons). The detection of y-ray s was soon recognised to be a 
promising channel (HQ]. If such a signal is yet to be measured 
by the existing y-rays observatories (FERMI, HESS, MAGIC, 
VERITAS), the prospect may significantly improve with the 
forthcoming next-generation instruments, such as CTA 0]. 

Having a modelling tool that can compute the expected y- 
ray flux from dark matter annihilation/decay in a wide range 
of astrophysical configurations and making it available to the 
community should prove very useful. This is the motivation 
for developing the code presented in this paper. This paper 
highlights clumpy's main features but a more thorough descrip- 
tion can be found in the documentation coming with the source 
code. The paper is organised as follows. Section [2] presents 
the formulation of the generic y-ray signal calculation that is 
performed by clumpy, with special emphasis on the contribu- 
tion from DM clumps. The code is described in sect. [3] which 
presents the structure, parameters and functions that are at the 
core of clumpy. Section [4] deals with two examples of the sci- 
ence that can be performed with clumpy. We conclude and dis- 
cuss the future developments of the code in sect. [5] 
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2. Calculating the y-ray flux from dark matter annihila- 
tion/decay 

The y-ray flux d<$> y /dE y from dark matter annihilat- 
ing/decaying particles is expressed as the product of a parti- 
cle physics term by an astrophysical contribution. For a given 
experiment with a spatial resolution a mt corresponding to an in- 
tegration solid angle AD = 2n(l - cos ai nt ), at energy E y and 
pointing in the direction {ip, 6), the flux is written as: 



squared along the line of sight (/, iff, ff), in the observational cone 



d<I> y 

-j^(£ 7 ,<M,AQ) = 



d<Pf 
dE v 



(E 7 ) X J(4r, 8, AD) . 



(1) 



2.1. Annihilation and decay: general statements 
2.1.1. Annihilation: the particle physics term 
It can be generically expressed as 



dd> 
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(2) 



where m x is the mass of the DM particle, <x ann is the annihila- 
tion cross section, (<x ann v) the annihilation rate averaged over 
the DM velocity distribution, Bf is the branching ratio into the 
final state / and dNl /dE 7 is the photon yield per annihilation. 
A widely accepted value for (cr ann v) estimated from present day 
thermal relic density is (cr ann v) = 3 ■ 10~ 26 cur's -1 [40. The 
photon spectrum is dependent on the annihilation channels. We 
restrict ourselves to the popular Minimal Supersymmetric Stan- 
dard Model (MSSM). In this framework, the neutralino is the 
lightest stable supersymmetric particle and is the favoured DM 
candidate. 

Neutralino annihilation can produce y-rays in three ways: i) 
DM annihilates directly into two photons or Z°y, giving rise 
to monochromatic lines, ii) or it can annihilate into primary 
products, the hadronization and decay of which produce a y-ray 
continuum, iii) finally, if charged particles are produced, inter- 
nal bremsstrahlung (IB) can also contribute to the continuum. 
The lines are very model-dependent and given the uncertain- 
ties on the particle physics models, we did not include any in 
clumpy. The user is however free to implement their own func- 
tion if need be. There exists in the literature several parametri- 
sations for the y-ray continuum and we have implemented that 
of Bergstrom et al. 0], and Tasitsiomi and Olinto The 
IB contributions is also very model dependent, and we provide 
one benchmark model from Bringmann et al. ffj. We refer the 
reader to 33.3.7l (spectra . h) for more details on these spectra. 

2.1.2. Annihilation: the astrophysical contribution 

The y-rays being produced from the annihilation of pairs of 
dark matter particles, the y-ray flux is proportional to the DM 
density squared. The second term of Eq. ([1), termed astrophys- 
ical contribution, corresponds to the integration of the density 



3 Note that a velocity dependent cross section should be included in a fu- 
ture release of the code in order to account for the Sommerfeld effect in the 
calculation. This effect results in much larger cross sections for some resonant 
neutralino masses and therefore boost the signal for these DM candidates. 



of solid angle AD (see Appendix A and geometry . h for defi- 
nitions and the geometry used) 



J(iff, 8, AD) 



Jl.o.s 



p 2 (l(if>,8)) d/dQ 



(3) 



where p(l(i//, ff)) is the local DM density at distance / from Earth 
in the direction iff, 8 (longitude and latitude in Galactic coordi- 
nates). 

This integral is the core of the clumpy code and can be rather 
complicated to evaluate. The difficulty comes partially from the 
steepness of some DM profiles considered but also, and mainly, 
from the existence of DM clumps, the contribution of which 
must be added to the smooth Galactic DM background. Both 
components are dealt with hereafter. 

2.1.3. Decaying dark matter 

y-rays can also be produced from the decay of DM. Formally, 
the flux is still given by Eq. ([]]) but the astrophysical and particle 
physics terms take different forms: 

• The particle physics term now corresponds to the decay 
spectrum and no parametrisation of that term has been 
included in the code so far. The parametrisation from 
Bertone et al. J3l could be used. 

• Because we are now considering a decay reaction, the as- 
trophysics term is the integration of the density (rather 
that density squared as in Eq. (01) along the line of sight, 
namely 



J (iff, 8, AD) 



Jo Jl.o 



p(l(if/,8)) dldD 



(4) 



This quantity can be straightforwardly com- 
puted in CLUMPY by setting the flag parameter 
gDM_IS_ANNIHIL_OR_DECAY to false (see table©. 

2.2. Dark matter distribution 

To date, there is no consensus to what the Galactic DM pro- 
file, p(r), should be. There is some dynamical evidence that 
the Galactic DM halo might be triaxial |@1 and numerical sim- 
ulations, such as the Aquarius |@] or the Via Lactea runs lioll 
also find non-spherical halos. When included, the baryons tend 
however to make the halos more spherical ll 111 . In this first ver- 
sion of clumpy we simplify the problem (as often done) using 
spherically symmetric DM halos, but allowing for tri-axiality is 
being considered for future releases. The total averaged den- 
sity profiles of the DM halos measured in these simulations are 
generally fitted by a form of the Zhao profile Il2ll : 



PtotW 



Ps 



(r/r s y [l+(r/r s )<*] 



a-\(p-y)l<* 



(5) 



where y, a and /3 are respectively the inner, transition and outer 
slope of the profile, p s and r s a scale density and radius. The 
Navarro et al. II 1311 . Moore et al. 11411 . and Diemand et al. 
| f5|] profiles fall in this category. Some other profiles having 
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a log-varying inner slope | If], 12 ] are a ls° included in clumpy. 
The reader is referred to 33.3.41 (profiles .h) for a detailed 
description. 

Some of these DM profiles are steep enough in the inner re- 
gions to lead to a singularity of the y-ray luminosity at the cen- 
tre of the halo. However, a cut-off radius r cut naturally appears, 
within which the DM density saturates due to the balance be- 
tween the annihilation rate and the gravitational infalling rate 
of DM particles. The saturation density reads 1 1811 



Psat = 3.10 1 



'( m * )x(l 

UOOGeV/ \ 



l(T 26 cm 3 s" 



<crv> 



M kpc" 3 . (6) 



Plugged in any profile parametrisation, this saturation density 
defines the cut-off radius below which the annihilation rate 
is constant. This is a very crude description, but this is not 
important as this cut-off matters only for the steepest profiles 
(y > 1 .5) which are disfavoured by current numerical simula- 
tions. 

The hierarchical scenario of galaxy formation in a A-CDM 
cosmology is characterised by a high degree of dumpiness of 
the DM distribution (supported by N-body simulations) so that 
the total averaged density computed from these simulations cor- 
responds to 



PtotO) = PsmOO + (PsubsM) 



(7) 



where p sm (r) is the "true" smooth component and (p su b(r)) rep- 
resents the average density from the substructures. These two 
components are discussed below. 

2.2.1. The clumps 

Whether substructures in subhalos are scaled-down versions 
of substructure in main halos remains an open question |@]. In 
this release of clumpy, we only consider one level of substruc- 
ture within the halo under scrutiny (Galactic halo or individual 
halo such as a dwarf spheroidal galaxy), and the properties of 
these sub-halos can be independently chosen from that of the 
parent halo (see below). 

Clump individual properties. In the code, all clumps have the 
same inner DM distribution, which can be any of the profiles 
discussed in 33.3.41 The mass of the clump generally suffices to 
determine all its properties, i.e., its size R~ k , and once an inner 
density profile is chosen, its scale radius r s and scale density 
p s . This is done using the so-called concentration parameter 
c v j r . Parametrisations of the mass-concentration relation have 
been established from numerical simulations for isolated field 
halos Il9ll20jl _but the concentration generally presents a signifi- 
cant scatter l2lll . These relations have been implemented in the 
code, regardless of that scatter or of the environment and for- 
mation history of the halos that can also affect the value of the 
concentration. We refer the reader to 33.3.61 (clumps .ti) and 
the Doxy gen documentation for more details. 

Spatial and mass distribution of the clumps. Following Lavalle 
et al. [220, spherical symmetry is assumed and the two distri- 



butions are assumed independent If Mot is the total number 
of clumps in a DM host halo, then the overall distribution of 
clumps is written as: 



d 2 N 
dVdM 



dP v (r) d>P M (M) 



dV 



dM 



(8) 



where the spatial and mass distribution are probabilities, re- 
spectively normalised as: 



dPyjr) 
dV 



dV = 1 



and 



dT M (M) 
dM 



dM = 1 .(9) 



Analysis of N-body simulations have shown the mass distribu- 
tion to vary as a simple power law 



ctPm(M) 
dM 



kM" % withff M « 1.9. 



(10) 



Given the mass and spatial distribution, the total number of 
clumps Mot can be determined in two ways. If the mass fraction 
/ of clumps and the mass of the host M lot halo are known, 



Mot = 



<M lcl > ' 



with (Mid) the average mass of one clump. An other approach 
is to know the number of clumps M) in a given mass interval 
[Mi , Mj\ and to calculate the total number of clumps as 



N 



Jmi 



dM 



Clumpy uses the first approach when dealing with substructures 
in individual host halos (that are not the Milky Way) and the 
second one for the Galactic halo clumps (numerical simulations 
have found about 100 clumps with masses above 10 8 M in 
Milky Way-like halos). clumpy allows to select any of the pro- 
files given in 33.3.4l (prof iles . h) for the clump spatial distri- 
bution, independently of the choice made for the total density 
profile. 

2.2.2. The smooth component 

From the spatial distribution dPyi^/dV one can define the 
average clump density (pa) used m Eq. ^} as: 



<PsubsW> = fM tot 



dPyjr) 
dV 



(11) 



For a given total density profile and clump spatial distribution, 
Eq. (|7l) then gives the smooth density profile as 



Psm(r) = PtotM - <Psubs(r)> • 



(12) 



Using this approach, rather than providing independently a 
smooth and a clump profile, ensures (by construction) that the 
total density profiles follows the one measured in N-body sim- 
ulations. This has been discussed in Pieri et al. M23 I. 



4 Keep in mind that such an assumption fails near the Galactic centre, where 
small clumps are expected to be disrupted by tidal forces. In this direction, the 
signal from the smooth contribution dominates anyway. 
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2.3. Calculating the annihilation J -factor 

The astrophysical contribution to the annihilation flux is for- 
mally given by Eq. ©. The J-factor is more explicitly written 

as: 



r An r'"'» if n ,1 

Jo Ju, 1 { i ) 



l l dldQ. . 



(13) 



where p sm is given by Eq. (TTZl i and the second term corresponds 
to the sum of the inner densities squared of all clumps ; con- 
tained with the volume element. The latter should not be con- 
fused with (psubsW)- Three terms arise from this equation: 



■/subs = 2j P c1 

"-'0 ^'min V | I 



dldCl 



(14) 



(15) 



and the cross product 



J cross— prod — 

2 1 /~^2X rfWn - (16) 



2.4. Statistical considerations 

Taking % = 1.9 as a canonical value for the slope of 
the mass distribution, with 100 clumps above 10 8 M Q as de- 
termined from N-body simulations, and a minimal theoretical 
mass of 10~ 6 M for the smallest ones, the total number of 
clumps in the Milky Way halo (i.e., for AQ = 4n) is estimated 
to be ~ 10 14 , which is computationally prohibitive. Even if we 
are interested in a 5° x 5° skymap only (i.e., AQ ~ 6 x 10" 3 ), 
the number of clumps is ~ 5 x 10 10 . This is too large a number 
to realise a statistical drawing according to the APy/AV distri- 
bution and then numerically integrate the J component for each 
of these clumps. However, their high number means that a sta- 
tistical average can be safely performed, as long as the variance 
is small. The clumpy code computes the y-ray flux from DM 
clumps by either randomly drawing clumps given user-defined 
spatial and mass distributions and inner DM profiles, or by cal- 
culating the mean contribution whenever the variance of this 
contribution is smaller than some user-defined value. 

2.4.1. Averaged quantities from the continuum limit 

Given the clump mass and spatial distributions given above, 
it is possible to define average values for the primary random 
variables of the problem (a clump mass and position) and the 
one deriving from them. The mean mass and the mean lumi- 
nosity are respectively defined by 



(M) 



U) 



f 

JM 

J 

Jm 



M - AV M 
M — -AM , 
AM 



AM 



(17) 



(18) 



where £(M) is the intrinsic luminosity of a clump given by 

£(M)= f (p cl ) 2 dV. (19) 

Jv„i 



The average of the distance to any power n is written as 

XAfi Wm„ AID 
J i"^i 2 didn. (20) 

Combining these, it can be shown that the average contribution 
from the clumps is expressed as 



(-/subs) - Ntt 



r r — *** f 

Jo J/ min dV J Mn 



~AM 



In the same spirit, in the continuum limit the cross-product be- 
comes 

(-/cross-prod) = 21 I p sm (p sahs )dldQ, . (22) 



The importance of the cross product can easily be identified 
when considering that the clump spatial distribution follows the 
total average density profile. In that case, (p su bs('")) = fptot(r) 
andp sm (r) = (1 - /)p to tW, leading directly to /cross-prod /-/tot = 
2/(1 - f). For a typical clump mass fraction / = 10%, 
this cross product amounts to 18% of the reference J-factor 
■Aot = f Pto t dld£l . This shows that when an averaged clump de- 
scription is used to integrate the signal along the line of sight, 
this cross product is not negligible and should be taken into 
account. In Appendix B it is however shown that on an in- 



dividual basis, the cross product of one clump in an underly- 
ing smooth Galactic component can be safely discarded. This 
will apply to the clumps that need to be statistical drawn in our 
model (see below). 



2.4.2. Variance of the clump contribution a 2 



clumps 



Along a given line of sight and for a given integration angle 
aim, the average flux between Z m j n and Z max given by Eq. (LTTT i 
can also be expressed as 



(/subs) = (NdKJia) , 



(23) 



where (N c \) is the average number of clumps in this volume 
(given the distance interval and mass range considered), and 
(J\ c i) the average flux of one clump (see also 12411 ). Using the 
point-like approximation, a clump of mass M at a distance / has 
a flux 



./ 



i a 



r- 



(24) 



so that, recalling that the mass of the clump and its location are 
independent variables, 



(25) 



The variance on the flux of a single clump is then defined as 

1 



The variance of a population of (N c \) clumps follows as 



<Af c i>cr? 



lcl ■ 



(26) 



(27) 




1 . (kpc) 

nun c 

Figure 1: Relative error between the average and true J value from clumps 
when the signal is integrated between / ra ; n and Z raax , according to Eq. (28), for 
M c \ e [1, 10] M Q . The signal is integrated over a mt = 0.1° (solid line) or 
a mt = 0.01° (dashed line). Three directions in the Galactic plane are consid- 
ered: towards the Galactic centre (blue), towards the anti-centre (black) and 
towards the Galactic East (red). For the latter the green dotted line illustrates 
how to get the critical distance (here 2 kpc) below which clumps (in this mass 
range) must be drawn if the user requires an accuracy of 5%. All distributions 
(smooth, clump distribution and inner DM clump profile) are taken as NFW 
profiles. 

2.4.3. Criterion for using the averaged description 

For a given integration domain, there is a threshold mass 
above which the clumps are not numerous enough to be de- 
scribed by an averaged description. Conversely, for a given 
mass decade there is a critical distance below which the aver- 
aged description fails and where clumps should be statistically 
drawn. 

The averaged description can be safely used, as long as the 
relative error RE made on the clump flux with respect to the 
total flux (i.e. including the smooth contribution) 

VWci <t 1c1 



RE j dum s = - 

c umps N cl (J 1 c l) + J smooth 

is smaller than some user-defined prescription. 



(28) 



Critical distance l cl i t for the clumps in the Galaxy. The RE is 
plotted in Fig. [T] as a function of the lower limit of the integra- 
tion / m ; n . As can be seen, for any required precision (e.g., 5% in 
Fig.Q] tagged by the dotted green line), one can infer the corre- 
sponding distance / cl j t below which the relative error becomes 
larger than the request. 

In practise, it is convenient to perform the calculation in 
terms of mass decades: for each decade, we find (by dichotomy) 
the critical distance below which clumps must be drawrjfl A 
lower limit of 10~ 3 kpc is set on l CI n to avoid divergences for 
a clump sitting at the observer's location. An example for the 
final number of clumps to draw below / cr ; t is provided in Tab. Q] 
The flowchart given in Fig.|2]summarises the various steps im- 
plemented in clumpy to obtain the skymap (see, e.g., Fig. [3}. 



Table 1 : Example of Z cr j, and the associate number of clumps to draw toward the 
anti-centre (criterion REj clum < 5%) for a 2VT skymap with an integration angle 
(instrument) of a mt =0.\° . Actually, a given / cl j t gives <n c i> number of clumps: 
na to draw is obtained from a Poisson distribution of mean value <n c \ >. 



Mass decade 
lo gio(t&) 



# clumps 
(Galaxy) 



Zcrit (kpc) 
kpc 



<# clumps> 

(5°x5°) 



# to draw 

(5° x5°) 



-6 


-5 


3.5 x 10 14 


1.0 x lO" 3 


1.5 x 10 -3 





-5 


-4 


4.4 x 10 13 


1.0 x 10" 3 


1.8 x 10- 4 





-4 


-3 


5.6 x 10 12 


1.0 x 10~ 3 


2.4 x 10~ 5 





-3 


-2 


7.0 x 10" 


7.3 x 10 -3 


1.2 x 10~ 3 





-2 


-1 


8.8 x 10 10 


5.8 x 10- 2 


7.6 x 10" 2 





-1 


: 


1.1 x 10 10 


4.0 x 10~' 


2.8 


3 





1 


1.4 x 10 9 


2.0 


3.7 x 10 +1 


45 


1 


2 


1.7 x 10 8 


6.7 


1.0 x 10 +2 


117 


2 


3 


2.2 x 10 7 


1.6 x 10 1 


8.7 x 10 +1 


99 


3 


4 


2.8 x 10" 


3.3 x 10 1 


3.6 x 10 +1 


44 


4 


5 


3.5 x 10 s 


5.8 x 10 1 


1.0 x 10 +1 


7 


5 


6 


4.4 x 10 4 


9.5 x 10 1 


2.2 


1 


6 


7 


5.6 x 10 3 


1.5 x 10 2 


4.1 x 10-' 





7 


8 


7.0 x 10 2 


2.5 x 10 2 


6.6 x 10~ 2 





8 


9 


8.8 x 10 1 


2.9 x 10 2 


9.5 x 10~ 3 





9 : 


10 


1.1 x 10 1 


2.9 x 10 2 


1.3 x 10~ 3 






This mass decade corresponds to a typical clump radius R c ] and we also 
check that the point-like hypothesis used to derived the continuum limit is ap- 
propriate (/ crit » R a ). 



Threshold mass m t hresh far the sub-clumps in an "far-away" 
halo. The second case where the criterion is used is for the 
sub-clumps in a halo located away from us, e.g., a dSph galaxy. 
In that case, the distance d to the clumps is known as they are 
all distributed within their parent halo's radius and we can de- 
termine in a similar fashion the threshold mass above which 
clumps must be drawn. 

3. Description of the code 

Clumpy is written in C/C++ and relies of the CERN ROOT 
library that is publicly available from h ttp : //root . cern . ch| 
The main features of the code are described hereafter along with 
brief descriptions of its most important functions. The code's 
structure is standard, with separate directories for the various 
pieces of code: declarations are in include/*. h, sources in 
src/* . cc, compiled libraries, objects and executables are re- 
spectively in the lib/, obj/, and bin/ directories. 

3.1. Code executables 

The code has one executable that can take three main op- 
tions depending whether the user i) is interested in performing 
a generic analysis on the Galactic halo, ii) is focusing on a list 
of specific halos (e.g. dSph galaxies), or iii) wishes to perform 
a statistical analysis on a single object. These main three usages 
of the code all accommodate several cases, the details of which 
are given below. 

3.1.1. Galactic mode: ./bin/clumpy -g [option] 

This calculates the /-factor for the Galactic halo, including 
both the smooth and the clump component. For the clump com- 
ponent, the mean flux or a stochastic realisation can be chosen. 
Note that a list of specific halos can be added to the calcula- 
tion. Several options allow to select the quantity to be calcu- 
lated/plotted from a text interface. A mandatory input is the 
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CLUMPS 



READ IN USER-DEFINED PARAMETERS 
DEFINE FIELD OF VIEW 



SMOOTH 



( COMPUTE L_CRIT(M) GIVEN CLUMP DISTRlF ) 



For a given mass decade, need to 
determine the distance l_crit below 
which the analytical clump estimation 
fails and clumps need to be drawn 
from dP/dV 



■a 

B 
■3 



COMPUTE NUMBER OF CLUMPS TO 
DRAW IN THE FIELD OF VIEW 
FROM to L_CRIT(M) 



COMPUTE ANALYTICAL CONTRIBUTION 
OF THE CLUMP "CONTINUUM" 



DRAW CLUMPS 



COMPUTE SMOOTH CONTRIBUTION 



IDENTIFY CLUMPS BELONGING 
TO THE F.O.V 



COMPUTE CLUMP FLUX AND 
STORE IN MAP PIXELS 



ADD ALL CONTRIBUTIONS PIXEL BY PIXEL 



STORE ALL COMPONENTS IN FILES 
AND/OR 
DISPLAY WITH ROOT 



Figure 2: Flow chart of the different steps leading to the creation of a y-ray skymap. 



Table 2: Format of the ASCII file to describe halos processed by . /bin/clumpy_dsph -h [option] . 'Type' should take one of the values of gENUM_TYPEHALOES, 
namely DSPH, GALAXY or CLUSTER, 'long', Tat' and 'd' are the galactic longitude, latitude, and distance to the halo, z its redshift and R v j, its virial radius. The profile 
parameters are p s , r s , 'prof.' (see prof ile .h) and the shape parameters. In particular, if 'prof'=kZHA0, see j)3.2K the profile corresponds to the (a, /}, y) profile 
(in which case the next-to-last three columns are used). If an kEINASTO is used instead, only the first shape parameter is considered. 



Name 


Type 


long. 


lat. 


d 


z 




Ps 


r s 


prof. 


a 


P 


1 






deg 


deg 


kpc 




kpc 


M G kpc~ 3 


kpc 











parameter file clumpy_params .txt described in §13. 21 and ta- 
ble |3] The options are the following: 

• 'option = -gl' plots the total, smooth and averaged 
clump density in ID, given the user-defined total density 
and the number of clumps in a given mass range. 

• 'option = -g2' plots J sm + <7 sub s> + <7 C ros S - P rodJ! re- 
spectively defined by Eq. ( fl4l i. Eq. (l2TT i and Eq. ( l22l as a 
function of the integration angle a mt . 

• 'option = -g3' is the same as '-g2' but as a function 
of the angular distance to the Galactic centre, for a given 
integration angle. 



Does not exist for decaying DM. 



• 'option = -g4' is the same as '-g3' but including a 
list of specific objects, such as dSph galaxies. 

• 'option = -g5' plots a 2D skymap of the /-factor, us- 
ing the averaged clump description, where 

• 'option = -g6' is similar ' -g5 ', but adds specific ha- 
los to the plots (e.g. dSph), to identify their contrast to the 
background. 

• 'option = -g7' and -g8 are respectively the same as 
' -g5 ' and ' -g6 ' but with the clumps randomly drawn 
from their user-defined mass and spatial distributions. 

3.1.2. Halo list mode: ./bin/clumpy -h [option] 

This usage of clumpy performs operations that are basically 
the identical to those performed in the 'Galactic' mode (see 
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{33.1 .It . but for a list of halos that are not the Galactic halo. 
This mode was initially implemented for the study of dwarf 
spheroidal galaxies 12511 . but can be extended to any DM halo 
(e.g., galaxy clusters). If only the smooth total density profile 
is considered, the halo profiles are read from a simple ASCII 
list of profiles along with their parameters, as exemplified in 
Tab. |2] If sub-clumps are considered, the clumpy parameter file 
clumpy_params . txt must also be provided (see table[3j. 

Again, further options are accessed for this mode by means 
of a simple text-user interface: 

• ' -hi ' to ' -h3 ' are identical to ' -gl ' to ' -g3 ' but trans- 
posed to any halo that is not the Galactic DM halo; 

• ' -h4 ' to ' -h.5 ' create 2D skymaps of the halo under 
scrutiny using respectively an averaged description or a 
statistical realisation for the sub-clumps within this halo; 

3.1.3. 'Statistical' mode: ./bin/clumpy -s [option] 

Confidence levels on the astrophysical factor J, for e.g. 
dSphs, can also be estimated from an ASCII statistical file. The 
format for the latter is quite similar to the halo list file (see doc- 
umentation) as it requires the profile parameters, but it differs 
in the sense that for a single halo, thousands of profiles 'se- 
lected' by a statistical analysis must be provided (their respec- 
tive x 1 value is necessary, although not used for all CL calcu- 
lations). Such a statistical analysis can be the Markov Chain 
Monte Carlo (MCMC), as discussed in Charbonnier et al. 112511 : 
Walker et al. 112611 . The most important options in the statistical 
mode are: 

• ' -s2 ' plots the probability density functions (PDF) and 
correlations of the parameters given in the statistical file; 

• ' -s4' locates the x% confidence intervals on the PDF of 
each parameters; 

• ' -s5 ' computes, for each radius from the dSph centre, the 
median and confidence levels of the DM density p{r); 

• ' -s6 ' is the same as ' -s5 ' but for J(a m ), i.e. as a func- 
tion of the integration angle; 

• ' -s7 ' is same as ' -s5 ' but for J(6), i.e. as a function of 
the angle to the centre of the halo. 

3.2. User-defined parameters 

For most of the options used above, a user-define parameter 
file is read. Table[3]gives the list of the main parameters. Each 
parameter is a single word formatted as g (for global parame- 
ter) + physics keyword + parameter name. For instance, COSMO 
refers to cosmology-related parameters, TYPE defines the type 
of halo under scrutiny (dSPh galaxy, galaxy or galaxy clusters) 
and GAL refers to the Milky Way halo parameters. Note also that 
the parameters ending in FLAG_CVIRMVIR, FLAG_SPECTRUM, 
and FLAG_PROFILE should take values from the following enu- 
merations: 

• enum gENUM.CVIRMVIR {kBOl, kENSOl, kJSOO) for 
the clump concentration; 



• enum gENUM.GAMMASPECT {kSUSY_BUB98 , 
kSUSY_T002 , kSUSY_BBE08) for the y-ray spectrum; 

• enum gENUM -PROFILE {kZHAO , kEINASTO , 
kEINASTOJJ} for the DM spatial distributions (smooth 
and clumpy components). The corresponding shape 
parameters (e.g. (a,/3,y) for kZHAO, a for kEINASTO and 
n for kEINASTOJJ) should be provided by the user in the 
input parameter file (see below). 

This enumerations allow the user to use their explicit denomina- 
tion kXXX, rather than the corresponding integer number. These 
parameters are gathered in the ASCII file clumpy_params . txt 
which is formatted as "Parameter name", "Units", "Value". 
The function loacLparameters (f ilejiame) defined in 
params.h reads in the input file and assigns the appropriate 
values to the parameters. The user can change any "Value" in 
this fileO. 

3.3. clumpy libraries 

The functions of clumpy are located in eleven "library" files, 
accordingly to their field of action. We briefly describe here- 
after the most important functions of these libraries. We are 
omitting in this document the more minor functions, the de- 
scription of which can be found in the documentation provided 
with the code. Also, we do not detail here what arguments are 
passed in these functions as it can be easily retrieved from the 
documentation. 

3.3.1. integr.h 

This library contains standard integration algorithms. The 
Simpson integrator with doubling step is the most widely used 
throughout the code to ensure that the user-defined precision is 
reached. Depending on the function to integrate, either linear 
or logarithmic stepping can be used. 

3.3.2. geometry. h 

In geometry. h are gathered functions performing geomet- 
rical tasks (change of coordinates, etc.). They are not essential 
for the presentation of the code and generally self-explanatory. 
We refer the reader to the extensive Doxygen documentation 



and to Appendix A for more details. 



3.3.3. inlines .h and misc .h 

These contain some definitions (unit conversion and inline 
functions) and miscellaneous formatting functions respectively. 

3.3.4. prof iles. h 

This library contains all the density profiles available in the 
code and operations on those. Among the functions in the "pro- 
files" section of the code, some act in a very straightforward 
way (e.g. squaring the density) and are not included in the 
present document. The reader is, once again, referred to the 
Doxygen documentation for a full description. 



7 Some of these parameters are optional and used only for specific 
modes [options] 
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Table 3: clumpy user-defined input parameter file (clumpy_params . txt) . 



Name 



Definition 



gC0SM0_RH00_C 

gC0SM0_0MEGA0_M 

gC0SM0_0MEGA0_LAMBDA 

gDM_GAMMARAY_FLAG_SPECTRUM 
gDM_MMIN_SUBS 
gDM_MMAXFACTQR_SUBS 
gDM_RHOSAT 

gTYPE_CLUMPS_FLAG_CVIRMVIR 1 
gTYPE_CLUMPS_FLAG_PROFILE 
gTYPE_CLUMPS_SHAPE_P ARAMS [0-2] 
gTYPE_DPDM_SLOPE 
gTYPE_DPDV_FLAG_PROFILE 
gTYPE_DPDV_RSCALE 
gTYPE_DPDV_SHAPE_P ARAMS [0-2] 
gTYPE_SUBS_MASSFRACTION 

gGAL_CLUMPS_FLAG_CVIRMVIR 

gGAL_CLUMPS_FLAG_PROFILE 

gGAL_CLUMPS_SHAPE_P ARAMS [0-2] 

gGAL_DPDM_SLOPE 

gGAL_DPDV_FLAG_PROFILE 

gGAL_DPDV_RS 

gGAL_DPDV_SHAPE_P ARAMS [0-2] 

gGAL_SUBS_Ml 

gGAL_SUBS_M2 

gGAL_SUBS_N_INMlM2 

gGAL_RH0S0L 

gGAL_RS0L 

gGAL_RVIR 

gGAL_TOT_FLAG_PROFILE 
gGAL_TOT_RSCALE 
gGAL_TOT_SHAPE_P ARAMS [0-2] 

gSIMU_ALPHAINT_DEG 
gSIMU_EPS 

gSIMU_IS_ANNIHIL_OR_DECAY 
gSIMU_SEED 



Critical density of the universe [M Q kpc ] 
Present-day pressure-less matter density 
Present-day dark energy density 

y-ray spectrum flag [gENUM.GAMMASPECT] 
Min. mass of a DM (sub-)clump [M G ] 
Defines max. mass of clump in host halo as M„ 
Saturation density [M Q kpc~ 3 ] 



Factor x M h , 



Clump concentration flag in the parent TYPE halo[gENUM_CVIRMVIR] 

Clump inner profile flag [gENUM_PR0FILE] 

Shape parameters for sub-clumps inner profile 

Slope subclump mass function 

Spatial distribution subclumps [gENUM_PRQFILE] 

Scale radius for the subclump distribution in the parent halo [kpc] 

Shape parameters for the clump spatial distribution 

Mass fraction of the parent halo in sub-clumps 

Concentration flag [gENUM.CVIRMVIR] 

Clump inner profile flag [gENUM_PR0FILE] 

Shape parameter for the Galactic clumps inner profile 

Slope the clump mass function 

Clump number distribution flag 

Scale radius for clumps [kpc] 

Three shape parameters for the Galactic clump spatial distribution 

Reference mass for gGAL_SUBS_N_INMlM2 [M G ] 

Reference mass for gGAL_SUBS_N_INMlM2 [M Q ] 

# of clumps in [M, , M 2 ] 

Local DM density [GeV cm" 5 ] 

Distance Sun - Galactic centre [kpc] 

Virial radius of the Galaxy [kpc] 

Total DM profile of the Milky Way [gENUM_PR0FILE] 

Scale radius for DM halo of the Milky Way [kpc] 

Shape parameter for the Galactic total density profile 

Integration angle [deg] 

Default precision for numerical integrations 

For annihilating or decaying DM 

Seed of random number generator 



TYPE corresponds to the type of halo under scrutiny (not the Galactic halo): DSPH, GALAXY or CLUSTER are the values 

allowed by gENUtLTYPEHALOES in this version. 
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Table 4: Standard DM density profiles parameters, following Eq. i29\ . 





a 


P 


r 


Reference 


ISO 


2 


2 







NFW97 


1 


3 


1 


Navarro et al. f 1 31 


M98 


1.5 


3 


1.5 


Moore et al. ri41 


DMS04 


1 


3 


1.2 


Diemand et al. 1151 



rho_ZHA0 returns the DM density of a user-defined profile 
having the generic functional form given in Zhao IU2I1 



P(r) 



ZHAO 



(r/r s y [\+(rlr s y] 



(29) 



The isothermal, Navarro et al. 1 131. Moore et al. lfl4ll . and 
Diemand et al. 115I1 profiles are all from the Zhao family. 
They are defined using Eq. (1291 with the set of parameters 
of table [4] The user is required to specify these parameters 
in the ASCII input file. 

rho_EINASTO and rho_EINASTO_N: it has been argued that 
a profile with a logarithmic inner slope gives a better fit to 
the simulated data. This so-called Einasto is generically 
written as 

p(r) cc exp (-Ar a ) . 



In Navarro et al. 1 16] and Springel et al. |s| the Einasto 
profile is defined as 



pEINASTO (r)= ( _- 

a 



— I -1 

r-2 



(30) 



where r_2 is the radius at which the profile's slope 
d In p/d In r = -2, and p_2 the corresponding density. 
Both Navarro et al. 13 and Springel et al. H find 
a ~ 0.17 to be a good fit to the smooth component of sim- 
ulated halos. Springel et al. Q further notes that a ~ 0.68 
allows a good description of the spatial distribution of sub- 
structures dP/dV. 

Merritt et al. Il7ll use an Einasto r 1 ^" profile, slightly dif- 
ferent from the above and given by 



p ElNASTO-N (r)=peexp 









-d n - 










) 



(31) 



where n = 6 and d„ = 3n — 1/3 + 0.0079/n. The scale 
radius r e in this relation correspond to the radius within 
which half of the halo mass is enclosed. 

The normalisation density of all these profiles (p s , p_2 
or p e ) can be done by i) either requesting a given value 
for the density at a given point or ii) by requesting 
a given mass within a given radius. The functions 
set_parO_givenjrttoref and set_parO_given_mref 
achieve the normalisation using both approaches. 



• rho simply returns the value of the density by calling any 
of the functions above, depending on the value of the pro- 
file flag given in clumpy _params . txt. 

• For two density profiles p\ and p2, rho_mix returns either 
the difference p\ - p2 or the product p\p2- This is needed 
when computing Eq. (fl2l i for the smooth density profile or 
Eq. ( l22l for the cross product. 

3.3.5. integr_los.h 

This is the core of the clumpy code as it contains the rou- 
tines performing the integration along the line of sight direction 
(0, if/) and over the chosen angular resolution a m . Given a func- 
tion g(l, a, j3\ 0, if/) where (I, a,/3) are spherical coordinates in the 
(0, if/) direction (see Fig. [6), the integration along the observa- 
tion cone, from l\ to h is generically written as: 

d/3 I sin a da I g(l, a,/3; 6, if/) fdl .(32) 
o Jo J/, 

• integrand_L: returns the value of the function 
/(/, a, (3; 0, if/) = g(l, a,/3; 0, if/) I 2 . A switch is included to 
integrate different functions of the density profile p. Note 
that here p can refer to any profile or distribution, be it the 
density or the spatial distribution of clumps dPy/dV de- 
pending on the normalisation used. The functions that can 
be integrated along the line of sight are: 

- /(/, a,/3; 0, if/) — p when calculating the clumps aver- 
aged contribution <7 su bs) of Eq. ( 1211 ). 

- /(/, a,/3; 0, iff) = p 2 in order to evaluate J such as in 
Eq. (Hl|. 

- f(l,a,/3;0,ifr) = l 2 p is used to calculate the mass 
fraction or number fraction of clumps in the integra- 
tion volume; 

- f(l,a,/3;0,if/) = Pp is needed when estimating the 
mean distance of the clumps as in Eq. ( 1201 ). 

- f(l,a,/3;0,iff) = l 4 p is needed when estimating the 
variance of distance of the clumps. 

- f(l,a,/3;0,ifr) = V 2 p must be evaluated to evaluate 
the variance on J as in Eq. 



Note that if decaying rather than annihilating DM is be- 
ing considered, i.e. gDM_IS_ANNIHIL_OR_DECAY is false 
instead of true, any p 2 is replaced by p in the above func- 
tion^. 

A variation of this function, called integrand_Lrel, is 
used whenever the steepness of the profile requires it (see 
Doxygen documentation for more details). 

f n_beta_alpha: computes the third integrand in Eq. (1321 1. 
namely the integration over I. The functions to inte- 
grate roughly behave as power laws, with a varying slope, 



This makes some of these functions redundant for decay, but they are still 
required for annihilation. 
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and log-step Simpson integration schemes are used (see 
integr.h). Some profiles can be very steep and several 
tests and tricks are needed to ensure the numerical stabil- 
ity/optimisation of the integration. The reader is again re- 
ferred to the documentation for an explicit description of 
the integration strategy. 

• f n_beta: computes the integration over the angle a in 
Eq. ( l32l ). A simple linear Simpson integration scheme is 
used. 

• los_integral: it is the final part of the triple integration, 
i.e. the integration over the angle B. This is again per- 
formed quite easily with a linear Simpson scheme. This 
function is called whenever a single density profile is to be 
dealt with. 

• los_integral_mix: as seen previously, it may be nec- 
essary to integrate a difference or a product of two 
densities along the line of sight. A flag parameter 
switch_rho is introduced to tell rho_mix to return ei- 
ther pi — p2 (switchjrho = 0) or p\p2 (switch_rho = 
1) (see prof iles .h). This function then behaves identi- 
cally to los_integral and uses the same flags with re- 
spect to the quantity to integrate. 

3.3.6. clumps. h 

All the functions related to the clumps mass and spatial distri- 
butions are gathered in this file. Here, we only give an overview 
of the most important ones. 

• add_halo_in_map: when using the code to generate a 
skymap where clumps are randomly drawn from their 
mass and spatial distribution, this adds the J contribution 
of the drawn clump to the appropriate pixels. The clump 
is only added to a pixel if 7 c i(pixel) > e x / sm (pix ) where 
pix () is the pixel containing the centre of the clump and 
e is the requested integration precision gSIMU_EPS. This 
avoids the clump from being integrated "too far" which 
would only be time consuming without changing the re- 
sult. 

• dpdv_XXX and set_normprob_dpdv: the former function 
returns the spatial distribution of the clumps as described 
in 32.2.11 and the second calculates the appropriate nor- 
malisation to make it a probability according to Eq. (0. 
The function dpdvJCXX is written in several versions de- 
pending of the reference frame and coordinate system in 
which the clumps are drawn (cartesian XXX=xyz, spherical 
centred on the galactic centre XXX=rthetapsi, spherical 
centred on Earth XXX=lthetapsi, or centred on Earth, in 
the cone AQ, in the direction (i/r, 6) XXX=lalphabeta). 

• dpdm and set_normprob_dpdm: this is the equivalent of 
the above but for the mass distribution of the clumps as 
given by Eq. ( TTUb . 

• f rac_nsubs_in_f oi returns the number fraction of 
clumps in a given field of integration w.r.t. the total num- 
ber of clumps in the host halo 



• f rac_nsubs_in_mlm2 returns the fraction of clumps with 
a mass in the range [m\,m2], given the total number of 
clumps. 

• j crossprod_contimmm returns the value of the cross- 
product integrated along the l.o.s as defined by Eq. (l22l . 
This is a simple call to los_integral_mix. 

• jcrossprod_interpolation interpolates the value of J 
in the pixels of a skymap, provided that / has previously 
been estimated in some locations of the map. This makes 
use of ROOT interpolation function based on Delaunay's 
triangles and allows to speed-up the calculation signifi- 
cantly when the skymap has a (too) large number of pixels. 

• j .smooth returns the /-factor when a single density 
profile is used (e.g. p tot ). This is a simple call to 
los_integral. 

• j_smooth_mix is the same as above but deals with a 
combination of densities, e.g. when computing Eq. (fT4l 
where p sm = p tot - (p c i). This is a simple call to 
los_integral_mix. 

• j-sub .continuum computes and returns the averaged J 
generated by the sub-clumps in a halo, in a given line of 
sight. 

• lum_singlehalo returns the total intrinsic luminosity of 
a DM halo, given its inner density profile as in Eq. ( fT9l ). 

• mass_singlehalo returns the total mass of a DM halo, 
given its inner density profile. 

3.3.7. spectra.h 

This library contains the different parametrisations of the av- 
erage y-ray spectra, 

dN y ^ dN f 

— - = > B f (33) 

dE y ^fdEy 1 

in the particle physics term of Eq. ©. As mentioned in 32.1.11 
we have included the Bergstrom et al.J3l and Tasitsiomi & 
Olinto parametrisations for the y-ray continuum and the 
Bringmann et al. |0] expression of the internal bremsstrahlung. 
Defining x = E y /m x , the three parametrisations are given by 
their corresponding functions below. 

• dNdE_BERGSTRDM returns 

dAV _ Q ?3 1 exp (-7.8jc) 
dE y m x x 1 - 5 

• dNdE_TASITSIOMI computes the y-ray spectrum as 

^r = J_ /!» 5 . 5 _5 -o.5 5 _ u \ 
dEy m x \ 3 4 2 12" / ' 
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dNdE_BRINGMANN estimates the internal bremsstrahlung 
contribution with 



dE v 



1 a 4(x 2 -x+l) 2 



n x(l — x + e/2) 



In 2 



l-x + e/2 



where a — 1/137 is the fine structure constant and e is the 
ratio of the mass of the W boson to the neutralino mass 
e = mw/m x . 

3.3.8. stat.h 

The functions included in this library perform simple statis- 
tical operations on the random variables of the problem (mass 
and distance of the clumps) and on quantities deriving from 
them (luminosity, J-factor). Ultimately, these are mainly used 
to determine the critical distance (or threshold mass) below 
(above) which clumps need to be drawn (see 32.4. 31 for the 
adopted procedure). Given these dark matter spatial (dPy/dV) 
and mass (df^M/dM) distributions defined in 32.2. Il l: 

• for a given mass and/or distance range, meanlcljnass, 
meanlcl, meanlcl_lum2, meanlcl_l, and meanlcl_j 
return respectively the average mass, luminosity, luminos- 
ity squared, distance and J-factor of one clump; 

• var lcljnass, var lcl_lum, varlcl_l, and varlcl_j re- 
turns the variance of the corresponding random variable; 

• finally, knowing the variance of the J-factor of one clump, 
f ind_lcrit_los and f ind_mthresli_los find the criti- 
cal distance and threshold mass for drawing the clumps. 
They both rely on a simple dichotomy to bracket the de- 
sired quantity, as illustrated on FigQ] 

3.3.9. j analysis. h 

This is the top-level library. It contains the main functions 
called when invoking any of the three options of clumpy (galac- 
tic halo, halo list or statistics, see 33.11 ). Note that the first word 
'xxx' of the function's name refers to which option they belong 
to (xxx=gal, xxx=halo, or xxx=stat). 

• xxx_j ID handles all the ID calculations of the Galactic or 
halo options. It computes p{r), J{ai nt ), and J(6), where 
is the angle measured w.r.t. the Galactic centre. 

• xxx_j 2D generates the 2D skymaps for the smooth, aver- 
aged and statistically drawn clump contributions. It also 
include the cross product term and It applies to both the 
Galactic and halo options. 

• xxx_load_list loads the list of halo under scrutiny either 
for the halo mode or the statistical option. The list should 
be provided as an ASCII file with the format given in ta- 
ble 

• xxx_set_par* are a set of functions filling the parame- 
ters arrays to be passed to the above functions with the 
proper quantities (e.g., profile flag and shape parameters, 
concentration flag, etc.). It applies to the Galactic and halo 
options. 



• stat_CLs is the most important function in the statistics 
option and computes the confidence levels of a quantity 
(density, J-factor) from its PDF. 

3 1 

x + x — — . 

4. Run examples 

This section provides a few plots obtained by using the 
clumpy executable. Many more examples and illustrations are 
provided along with the Doxygen documentation. 

4.1. Skymaps 

Figure [3] shows two examples of a 5° x 5° skymap of the 
J-factor (option ' -g7' of clumpy_gal), when looking 10° off 
the Galactic centre (left) or towards the anti-centre (right). The 
angular resolution is 0.01° and the Galactic smooth DM density, 
clump distribution and clump inner profile are all assumed to 
follow an NFW parametrisation. 

Similarly, skymaps of individual objects can be performed 
(option ' -h5 ' of clumpy_dsph). In Fig. [4] a "generic" dwarf 
galaxy, including a population of subclumps, is integrated over 
0.1° (left), 0.05° (middle), and 0.01° (right). As the resolution 
increases, more substructures become apparent but the J-factor 
is strongly reduced. 

4.2. dSphs and statistical tools 

The three panels in Fig[5]illustrate clumpy capabilities on the 
dSphs. The top left and bottom figures are taken from Char- 
bonnier et al. 12511 . and the top right one is a direct result of the 
analysis performed in the above paper. Running clumpy on the 
default parameters will not give these plots, but similar ones in 
spirit. 

5. Conclusions 

The clumpy code is optimised to calculate the line-of-sight 
integration of the density and/or the density square of the DM 
distributions (smooth and clump distributions). It can be used 
by (i) experimentalists looking for realistic y-ray skymaps to 
calculate their new instrument sensitivity, or simply to use them 
in model/template analyses for the DM diffuse emission near 
the Galactic centre (or for dSphs), (ii) for astrophysicists work- 
ing on the reconstruction of the DM content of dSphs who wish 
to calculate the J factor as a by-product of their study, (iii) 
theoreticians who want to plug their preferred particle physics 
model and see what is the corresponding y-ray flux in the 
Galaxy/dSph. The code is fully documented and provides many 
functions that makes it versatile and flexible. 

Many improvements are possible for a future release. The 
simplest and more straightforward ones to consider are: 

• Extra-galactic signal from the local structures (e.g. local 
clusters, Nezri et al. in prep.), but also from the high- 
redshift isotropic contribution. 

• The possibility to have the so-called Sommerfeld-like en- 
hancements. This requires to couple the calculation of the 
spatial density of DM to the position-dependant velocity 
of the DM particles. 
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Figure 3: Example of two 5° X 5 skymaps, at 10° from the Galactic centre (left) and towards the anti-centre(right), for a 0.01° angular resolution, obtained from the 
skymap mode of clumpy. The colour scale gives the J values. These maps are obtained with the following command option: . /bin/ clumpy -g7. Clumpy will 
display the results automatically but these images have been produced with Gnuplot for aesthetic reasons. 





I 



Figure 4: Skymaps of a fake dSph galaxy at three different resolutions: 0.1° (left), 0.05° (middle), and 0.01 "(right). Note that ROOT will automatically display the 
image, but these ones have been plotted using the output files with Gnuplot. These maps are obtained with the following command option: . /bin/clumpy -h5 



• Neutrino annihilation spectra and fluxes. 

Some other improvements are intrinsically more difficult to add, 
as they depend on astrophysical inputs that are not well deter- 
mined. Indeed, in the absence of a smoking gun signature for 
indirect detection (e.g. a line detection), multi-wavelength and 
multi-messenger analyses seem to be a compulsory approach. 
To do so, one needs to take into account 

• the anti-proton and anti-deuteron signal; 

• electron and positron production; 

• secondary multi-wavelength emissions from leptons. 

Finally, the clumpy code could be coupled to some parti- 
cle physics code to get more specific annihilation spectra (e.g., 
darkSUSY, Gondolo et al. Jzl; MicrOMEGAs, Belanger et 
al. 128[), or to some cosmic-ray propagation code (e.g. GAL- 
PROflE DRAGON0 or USINE, Maurin et al., in preparation) 



http://galprop . Stanford, edu/ 



http://www.desy.de/ maccione/DRAGON/ 



to compare with the Galactic astrophysical backgrounds. This 
will be considered for a second release of the code. 
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Appendix A. Geometry and change of coordinates 

We define two frames of reference: one attached to the ob- 
server the origin of which is the Earth, and the other linked to 
the Galactic Centre (see Fig. [6}. The frame related to the Earth 
(E, 1, if/, 6) is spherical, and is identified to the Galactic coordi- 
nate system. / is the distance to the Earth, the Galactic longitude 
if/ is measured in the plane of the Galaxy using an axis pointing 
from the Earth to the Galactic Centre, the Galactic latitude is 
measured from the plane of the Galaxy to the object. 
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Figure 5: Top left: boost factor as a function of the integration angle for 'several' generic dSphs ( . /bin/clumpy -h.2). Top right: median (solid line) and 95% 
confidence levels (dashed lines) of the /-factor as a function of the integration angle for the Draco dSph ( . /bin/ clumpy -s6). Bottom panel: the ./-factor for the 
dSph can be over-plotted on the galactic smooth and mean-clump / factors ( . /bin/ clumpy -g4). 



Line of sight (l.o.s.) direction (Y,0 ) 



Galactic Centre frame (Oxyz) 

- GC: (0, 0, 0) 

- Earth: (-Rsol, 0, 0) 

Galactic Coord. (1, V, 6 ) 

- Earth: (0, 0, 0) 

- GC: (Rsol, 0, 0) 




-W0 



Figure 6: The geometry as defined in the clumpy simulation. Two frames are used: the Galactic Centre frame and the Galactic coordinates. For a given direction of 
observation (if/, 0) (l.o.s. for line of sight), we will run over all the a and /? directions within the solid angle An. 
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Figure 7: Relative importance of the cross-product w.r.t to the total J calculation for i) a clump embedded in the Galactic halo in the anti-centre direction (top panels) 
and ii) a sub-clump in DM halo located 100 kpc away (bottom panels). Several masses for the clump have been considered and the integration angle is a-i llt = 0.1°. 



The Galactic Centre frame (GC, x, y, z) is a cartesian frame, 
where the (GC, x) and (GC, y) axes are in the Galactic plane, 
and the (GC, z) axis is pointing to the North Galactic pole. In 
this frame, the position of the Earth is (-R©, 0, 0), where R is 
the distance between the Earth and the Galactic Centre. 

The change of coordinates formulae are given in the docu- 
mentation (geometry . h). 

Appendix B. Cross-product of individual clumps 

In §12.3l it has been shown that the cross-product between the 
underlying smooth DM density and that of all the clumps in 
the line of sight (when using an average description) can sig- 
nificantly contribute to the value of J in the line of sight. It is 
interesting to check what contribution it brings when consider- 
ing one clump only (rather that the entire distribution). We use 
the same equations as in 92.31 but written for one clump only 
and numerically integrates J. Two configurations are explored: 
i) a clump embedded in the Galactic halo (fig. [7] top panels) 
and ii) a clump embedded in a host halo located away from the 
sun (fig. El bottom panels). 

In fig. |7] the relative error made when neglecting the cross- 
product is plotted against i) the distance to the clump for the first 



configuration (top), ii) the distance of the clump to the centre of 
its host halo (bottom). In both cases, these curves show that, as 
far as one individual clump is concerned, the cross-product is 
completely negligible, whatever situation occurs: it amounts at 
most to ~ 1%. 
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